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Abstract 

We describe a way to introduce physics high school students with no background in programming to computa- 
tional problem-solving experiences. Our approach builds on the great strides made by the Modeling Instruction 
reform curriculum. This approach emphasizes the practices of "Developing and using models" and "Computational 
thinking" highlighted by the NRC K-12 science standards framework. We taught 9th-grade students in a Modeling- 
Instruction-based physics course to construct computational models using the VPython programming environment. 
Numerical computation within the Modeling Instruction curriculum provides coherence among the curriculum's dif- 
ferent force and motion models, links the various representations which the curriculum employs, and extends the 
curriculum to include real-world problems that are inaccessible to a purely analytic approach. 



1 Introduction 

Numerical computation (the use of a computer to solve, 
simulate, or visualize a physical problem) has fundamen- 
tally changed the way scientific research is done. Systems 
that are too difficult to solve in closed form are probed 
using computation (e.g., models of cardiac arrhythmias 
[1]); experiments that are impossible to perform in the 
laboratory are studied numerically (e.g., the dynamics of 
black hole collisions [2]). Consequently, in modern science 
and engineering, computation is widely considered to be 
as important as theory and experiment. 

Unfortunately, most high school students today are 
never introduced to computation's problem-solving pow- 
ers. Computer usage is widespread in high school STEM 
courses (e.g., obtaining lab data using computer acquisi- 
tion hardware/software), but such usage rarely involves 
students constructing a computational representation of 
a STEM problem. The lack of computation in domain- 
specific STEM courses is not addressed in most high 
school computer science courses, which typically focus 
on programming and procedural abstractions rather than 
solving science problems. In recognition of these short- 
comings, the recently published National Research Coun- 
cil's (NRC) framework for next-generation K-12 science 
standards lists "computational thinking" as one of the 
fundamental "practices" that should be incorporated into 
future K-12 science curricula [3]. The framework acknowl- 
edges that experience with computational thinking is cru- 



cially important, not only for developing future scientists 
and engineers, but also for providing all citizens with gen- 
eral insight into the science behind proposed solutions to 
technically complex social problems. 

In this article, we describe a way to introduce physics 
high school students with no background in programming 
to computational problem-solving experiences. Our ap- 
proach builds on the Modeling Instruction reform cur- 
riculum, which is currently used in approximately 10% 
of U.S. high school physics classrooms [4]. The Mod- 
eling Instruction approach emphasizes the practice of 
"Developing and using models" highlighted by the NRC 
K-12 science standards framework [3]; coupling com- 
putational experiences with Modeling Instruction en- 
ables the model-development practice and the compu- 
tational thinking practice to reinforce each other. To 
achieve this synergy, we taught 9th-grade students to use 
the VPython programming environment [5, 6] within a 
Modeling-Instruction-based physics course [7]. We found 
that numerical computation within the Modeling Instruc- 
tion curriculum provides coherence among the different 
models within the curriculum, links the various repre- 
sentations that the curriculum employs, and extends the 
curriculum to include real-world problems that are oth- 
erwise inaccessible using a purely analytic approach. We 
briefly describe the Modeling Instruction curriculum for 
the benefit of the reader, highlighting the areas where 
numerical computation can be embedded. We describe a 
sample activity in which 9th grade students have engaged 
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for concreteness. Finally, we reflect on our experiences, 
emphasizing our successes and our challenges. 

2 Modeling Instruction & 
Numerical Computation 

The Modeling Instruction curriculum employs a coherent 
framework of scientifically testing the limits of physical 
models (i.e., "the modeling cycle") by engaging students 
in the construction and comparison of different represen- 
tations of physical phenomena. Because of its emphasis 
on models, its focus on inquiry, and its use of multiple 
representations, the Modeling Instruction curriculum is 
effective not only in teaching students physical concepts 
[8], but also in encouraging participation in class [9], in 
helping align students' views about the nature of science 
with expert views [10], and in promoting students' self- 
efficacy [11]. 

In the Modeling Instruction curriculum, students con- 
struct models through the modeling cycle. Each mod- 
eling cycle is built on a set of modules; these modules 
promote scientific thinking through observation, experi- 
mentation, and discourse. By observing physical phenom- 
ena, representing those phenomena in a variety of ways, 
and making predictions of similar but not-yet-observed 
phenomena, students construct a working model which 
is able to fully describe the phenomena they observe. In 
the Modeling Instruction curriculum, the word "model" 
is typically used to describe the simplified mathematical 
representation of a physical system. For example, after 
observing a cart moving on a track at constant velocity, 
students describe the motion of the cart with a story (a 
linguistic representation). Kinematic graphs (graphical 
representations) are introduced as a useful representation 
of the motion. Students then connect these graphs to a 
motion map, a diagrammatic representation that uses ar- 
rows to describe the velocity of the cart along the track. 
This model is then extended to other physical situations 
for which the model applies. Students use this model to 
make predictions and draw conclusions based on these 
multiple representations. The instructor then presents a 
new physical system for which the developed model is in- 
sufficient (e.g., an accelerating car cannot be described 
using a model for constant velocity). When a model can 
no longer describe the phenomena being observed, the 
model is extended and the cycle is repeated (e.g., by in- 
vestigating a model for constant acceleration). 

Modeling Instruction treats each force and motion 
model as distinct, but the common thread of predict- 
ing motion using Newton's 2nd law and kinematics uni- 
fies them. The computational algorithm used to pre- 
dict motion likewise retains the distinctions between the 
force and motion models, but highlights the commonality 
among them: namely, that such models differ only in the 
net force exerted on the system and in their particular 
initial conditions. Given knowledge of the system's ini- 



tial position and velocity, as well as the net force on the 
system, the algorithm for predicting motion can be de- 
scribed as a set of rules applied locally in space and time: 
(1) At a given instant in time t, compute the net force, 
Fnet, acting on the system, (2) For a short time At later, 
compute the new velocity of the system using Newton's 
2nd law, (3) At the same new time (t + At), compute 
the new position of the object using this updated veloc- 
ity, and (4) Repeat Steps (l)-(3) starting at the updated 
time t + At. Formally, the iterative application of Steps 
(l)-(3) is, in effect, explicit (Euler-Cromer) numerical in- 
tegration [12] of the equations of motion for Newtonian 
mechanics (Av = a = F net /m, Ax — vAt). 

The mathematics behind iteratively predicting motion 
in this manner is well within the capabilities of most 
high school physics students (in either algebra-based or 
calculus-based courses); arguably, it is more accessible 
mathematically to students than the analytic methods 
currently used, even for the simplest cases (e.g., constant 
acceleration motion). Iterative motion prediction is usu- 
ally too labor-intensive to perform by hand, but these 
calculations can easily be handled by a computer. More- 
over, this same computational algorithm (Euler-Cromcr 
integration) can be used simulate the vast majority of 
physical systems at a high-school level, further reducing 
the barrier for introductory students to explore complex 
systems. 

Numerical computation offers significant pedagogical 
advantages. Computation highlights the relationship 
between the different physics models in the Model- 
ing Instruction program (e.g., the no-forces model, the 
balanced- forces model, and the unbalanced-forces model) . 
To produce simulations with qualitatively different be- 
havior, we simply change the initial conditions (e.g., 
from ID to 2D motion) or the net force (i.e., from con- 
stant to constantly changing). For example, we can 
generalize the constant velocity/balanced forces model 
to the constant acceleration/unbalanced forces model by 
inserting a constant net force into the computational 
model. Furthermore, we can extend the constant acceler- 
ation/unbalanced forces model to parabolic motion model 
by giving the object and initial velocity in both x and y 
directions. 

Numerical computation provides dynamic animation 
and visualization of representations that are otherwise 
static in the Modeling Instruction curriculum. The out- 
put of numerical computation is continuously-updating 
graphs (analogous to a chart recorder) and animations, 
not just numbers. The visualization provided by a numer- 
ical model is of paramount importance; certain aspects of 
visualization help students communicate a more coherent 
picture of their understanding [13]. These graphical and 
diagrammatic descriptions of the physical model (which 
might otherwise form the sole basis of the students' expo- 
sure to the model) are reproduced precisely by the com- 
putational model. Furthermore, the linking of represen- 
tations can be done quite easily with a few simple lines 
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of code (see Section 3) . 

These numerical models are not limited to analytically 
tractable solutions. This allows students to explore their 
real-world, rather than laboratory-constructed, observa- 
tions. Numerical computation provides a platform to fo- 
cus class discussion on modeling and investigation with- 
out the undue burden of sophisticated mathematical tech- 
niques. For example, students observe objects that expe- 
rience drag in their daily lives (try kicking a soccer ball 
on a windy day!), and yet a model of this phenomenon 
is not explored in most introductory physics courses. A 
model of turbulent drag is a simple model to construct 
and describe. We have found that students can construct 
a model for drag, make sense of the model's predictions, 
and compare those predictions to those of the constant 
acceleration model (see Section 4). 

By learning to use numerical computation, students ac- 
quire a familiarity with the tools of modern scientists and 
engineers. Moreover, as students gain experience with nu- 
merical computation, they begin to build "computer mod- 
els" as part of their normal practice of constructing and 
testing models, which further emphasizes that models are 
what scientists and engineers use to describe physical phe- 
nomena. Numerical computation can be an effective tool 
for exploring the limits and refining the physical model in 
question. Students can explore the influence of such in- 
puts on the resulting motion by changing parameters in 
the computational model. They can engage in prediction 
and confirmation by reviewing the animation and graph- 
ical output that their computational model produces. 

3 Developing a Set of 
Computational Tools 

Numerical computation can provide additional benefits to 
student understanding of science and, in particular, phys- 
ical phenomena. However, the tools that students use to 
numerically model such phenomena must include no more 
programming than is necessary. Their physics class is not 
a computer science course; hence, the program statements 
that students write should only reflect the representations 
with which they are becoming familiar. This is not to say 
that learning to use a particular programming environ- 
ment is not a learning goal for the course, but that we 
limit the amount of programming each student has to 
learn in order to reduce the burden on them. For our 
implementation, we developed a lightweight but powerful 
environment and employed it to focus students' computa- 
tional model development on the physics of the particular 
system and the representations of that model. We chose 
to use the VPython programming environment [6] along 
with a module developed for enhancing aspects of the sim- 
ulation, PhysUtil [14]. PhysKit [15], a web-based solution 
for installing and updating VPython and PhysUtil, was 
recently developed by team of Georgia Tech computer sci- 



ence seniors. This software is publicly available 1 and is 
in the testing phase. 

VPython is based on the Python programming lan- 
guage and provides a environment to write simple pro- 
grams that yield robust three-dimensional simulations 
(Figure 1). It has been used in conjunction with the 
Matter and Interactions curriculum [16] at a number uni- 
versities, colleges, and high schools [17]. The VPython 
programming environment was designed to limit the pro- 
grammatic statements needed to generate highly visual 
three-dimensional simulations. Students who receive suf- 
ficient computational instruction using VPython are able 
to successfully model novel situations [18, 19]. 

Figure 2 shows sample VPython code which models the 
motion of a fan cart subject to a single constant force. To 
construct this model, 9th grade physics students created 
the objects and assigned their positions and sizes (lines 
6-7), identified and assigned the other given values and 
relevant initial conditions (lines 9-10 and 12-14), calcu- 
lated the net force acting on the object of interest (line 
23), and updated the velocity and position of this ob- 
ject in each time step (lines 24-26). This code illustrates 
the algorithm students are taught to predict the motion 
of objects given the model for their interactions (Euler- 
Cromer integration) [12]. The code shown in Figure 2 
produces a highly visual simulation generated from a few 
program statements. This program represents what stu- 
dents are able to construct after instruction in our 9th 
grade conceptual physics course. 

The program shown in Figure 2 makes use the PhysUtil 
module. Developed as a team of Georgia Tech computer 
science majors' final design project, the PhysUtil mod- 
ule [14] was designed to further limit the code-writing 
needed to create highly visual simulations and to en- 
hance the functionality of VPython to include features of 
the Modeling Instruction curriculum (e.g., motion maps). 
The construction of the module focused on the develop- 
ment of Python classes that included additional represen- 
tations (e.g., motion maps and graphs) without the addi- 
tional burden of writing complex program statements. At 
present, we have added four Python classes with PhysU- 
til: PhysAxis, PhysTimer, MotionMap, and PhysGraph. 
The first three classes enhance the diagrammatic repre- 
sentation of the model by providing a baseline for length 
(PhysAxis) and time (PhysTimer) measurements as well 
as a representation of physical quantities (MotionMap) 
(Figure 1). The last class (PhysGraph) provides a graph- 
ical representation of the motion; physical quantities can 
be plotted as a function of time. Each of these classes re- 
quires a single initialization line (lines 16-19 in Figure 2), 
which can be provided to the students, and a single up- 
date line in the calculation loop (lines 28, 31, and 33). 
PhysUtil is currently maintained by the Georgia Tech 
Physics Education Research group and is available for 
free. 2 

1 https : / /per . gatech. edu/physkit/ 
2 https : / /per . gatech. edu/physutil/ 
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Figure 1: The visual output of a VPython+PhysUtil model of a soccer ball kicked in the air (without drag) constructed 
by three 9th grade students. PhysTimer appears in the upper-right corner (blue text). PhysAxis appears under the 
ground (blue line and text). MotionMap generated the "breadcrumbs" for the motion with time stamps and integer 
ordering (red spheres and red text). 



To illustrate how our particular brand of numeri- 
cal computation fits into a typical Modeling Instruction 
course, we present an activity used in a 9th grade physics 
course during the second half of the semester. Students 
employed and extended the parabolic motion model for 
the motion of an Angry Bird 3 to characterize the motion 
of a kicked soccer ball. 

4 A Typical Activity: Modeling a 
kicked soccer ball 

In our modified Modeling Instruction course, we pre- 
sented projectile motion after students had studied 5 pre- 
vious models (constant velocity, balanced forces, constant 
acceleration, unbalanced forces, and conservation of mo- 
mentum). Students discovered, through observation, that 
the constant acceleration model was insufficient to de- 
scribe the motion of objects in two dimensions subject to 
the ordinary gravitational force, F grav = mg. In fact, an 
appropriate description required the use of two models: 
the constant acceleration model in the vertical direction 
and the constant velocity model in the horizontal direc- 
tion. Typically, the parabolic motion model represents 
the capstone of the Modeling Instruction curriculum's 
treatment of force and motion. Later in the course, stu- 
dents explored motion confined to a circle by applying the 

3 http : //www.rovio . com/ en/ our- work/game s/view/1/ 
angry-birds 



circular motion model, but the parabolic motion model 
was the only true multidimensional model presented in 
the curriculum. In our treatment, we used numerical 
computation to investigate the parabolic motion model, 
to compare its predictions to real-world observations, and 
to resolve the limited predictions of this model by extend- 
ing the model to include air-resistance drag. 

Typically, students collected data from a lab experi- 
ment to motivate the development of a new models in 
the Modeling Instruction curriculum. It is also possible 
to collect data from msomething which is itself a model 
(e.g., a computer game). We motivated the parabolic 
motion model by showing students a snapshot of the tra- 
jectory of an bird from the popular Angry Birds video 
game (Figure 3). From this vector construction, students 
concluded that there must be a force acting on the Angry 
Bird which points vertically downward. To investigate 
this claim, we collected video data of Angry Birds flying 
across the screen using QuickTime's screen capture capa- 
bilities. The video data was then imported into Tracker, a 
free and open-source video tracking software [20], where 
the motion of the Angry Bird was logged and plotted. 
Tracker allows the user to plot a variety of quantities (i.e., 
x versus y, x versus t, et cetera). It can also compute, 
using discrete derivatives, the velocity and acceleration 
of the tracked particle in each coordinate's direction and 
plot those quantities. It also possible to do this analysis 
using Python, but the current implementation requires 
far too much code for students to handle reasonably. We 
are working on solutions for this. Students concluded that 
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1 from future import division 

2 from visual import * 

3 from physutil import * 

4 from visual. graph import * 

6 track = box ( pos=vector ( , -0.05, 0), size=(5.0, 0.05, 0.10), 

color=color .white ) 

7 cart = box(pos=vector(-track. length/2, 0, 0), size=(0.1, 0.04, 0.06), 

color=color . green ) 

8 



12 deltat =0.01 

13 t = 

14 tf = 6.45 
15 

16 timerDisplay = PhysTimer ( 1 , 1) 

17 graph = PhysGraphf ) 

18 axis = PhysAxis (track, 10, axisColor=color . red) 

19 motionMap = MotionMap( cart , tf, 10, markerType= " breadcrumbs " , 

labelMarkerOf f set=vector ( , . 3 , ) , dropTime=True ) 

20 

21 while t < 2.3: 

22 

T3 Ffan = vector (-0.75, 0, 0) 

24 accel = Ffan/mcart NeWtOtl's Second Law 

25 vcart = vcart + accel * deltat + Kjnematics 

26 cart.pos = cart.pos + vcart * deltat 
27 

2 8 motionMap . update ( t ) 

29 

30 t = t + deltat 

3 1 timerDisplay . update ( t ) 
32 

33 graph. plot ( t , cart.pos.x) 

Figure 2: A student's VPython program that models 
the motion of a fan cart subject to a constant force 
(constant acceleration/unbalanced forces model). Green 
boxes highlight where we focus students' attention during 
model construction. 

the constant velocity model applied to the bird's horizon- 
tal motion and the constant acceleration model applied 
to the vertical motion. These two models connected in 
this way compose the parabolic motion model. 

In our course, students also generalized the constant- 
acceleration computational model which they had devel- 
oped to the parabolic motion model. Students had pre- 
viously developed a fully "vector compliant" program to 
model constant acceleration in either the x or y direction 
(e.g., lines 23-24 in Figure 2). This generalization em- 
phasizes the interconnected nature of different force and 
motion models in the Modeling Instruction curriculum. 
The generalization is quite simple because, in a compu- 
tational model, the change beween different types of mo- 
tion under constant acceleration is simply a change of the 
initial conditions. By giving the object an initial velocity 
with a nonzero component orthogonal to the acceleration, 
the student can very easily move from modeling an object 
dropped from a known height (for example) to an object 
fired into the air at a known angle. 

The new computational model linked the different rep- 
resentations of the physical system and provided instant 
visual feedback about students' physical model. Using 
PhysAxis, PhysTimer, and PhysGraph, students recon- 
structed the Angry Bird's motion. Moreover, a compu- 
tational model allowed students to immediately observe 
if their physical model had any inconsistencies (e.g., un- 
expected motion in the horizontal direction) or if their 
computational model had any unrealistic effects (e.g., mo- 
tion not terminating at the ground level). The latter led 
to a nice discussion of the limitations of computational 
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Figure 3: A vector construction to determine the direc- 
tion of the acceleration in the Angry Bird's world. Fair 
use reproduction (non-profit educational illustration). 

models; they can only do what you have told them to 
do. Using MotionMap, students constructed an animated 
motion map to observe how components of the force or 
velocity change with time. With a computational model, 
students were able to systematically adjust parameters 
(e.g., the Angry Bird's mass, size, and initial velocity) 
to observe their effects on the animation; students paid 
particular attention to their graphs of kinematic and dy- 
namic quantities and their motion maps. Students re- 
ported their observations to their peers. 

Students were then confronted with the following chal- 
lenge: "We have learned that the constant acceleration 
model can help us describe how an object moves in one 
dimension, and that the parabolic motion model can help 
us in two dimensions. Now, what about a soccer ball that 
you kick into the air? How can we model this situation?" 
Typically, this would be dealt with using the parabolic 
motion model. By using the computational modeling, we 
can push this further. "What about a real kick, outdoors? 
Maybe it's a windy day. Do any of these models still ap- 
ply to the motion?" 

Students concluded from video analysis of a kicked ball 
(similar to the Angry Birds analysis) that there were ac- 
celerations in both the horizontal and vertical directions. 
Moreover, they observed that these accelerations changed 
with time. Students proposed air resistance as the cul- 
prit for this change. However, the model for air resistance 
(even linear drag) does not lend itself to analytical solu- 
tions achievable by 9th grade students; the mathematics 
is too sophisticated. Computational modeling allowed us 
to insert a velocity dependent drag force on the ball, and 
then to use the model to predict the trajectory and land- 
ing point for the soccer ball. 

Students further generalized their computational model 
to investigate features of this motion using the tools of 
PhysUtil. Using MotionMap, students observed that the 
direction of the force changed throughout the motion. 
Moreover, they discovered that the velocity component 
in the x direction was always reduced as the ball moved. 
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Using PhysAxis, PhysTimer, and PhysGraph, they were 
able to reproduce the motion of the ball with reasonable 
agreement with their observations. The motion of real 
projectiles is no longer intractable to conceptual physics 
students. 

Using numerical computation in the way we described 
does not supplant the typical activities in which students 
engage; it enhances and extends those activities. We are 
not replacing the great strides made by the designers of 
the Modeling Instruction curriculum; rather, we are en- 
hancing their worth by making the activities more rele- 
vant to students by including real-world examples, em- 
phasizing the concept of models, illustrating the gener- 
ality of physical principles, and providing a platform for 
future learning in numerical computation. 

5 Reflections 

We have used numerical computation in a 9th-grade 
Modeling-Instruction-based Honors Physics course in a 
private school setting for the last two years, each year 
comprising a different set of 15-18 students. In that time, 
we have observed several challenges to student learning 
and broader adoption. 

Students find debugging their programs difficult; that 
is, they have trouble determining whether they have made 
a coding error or a physics error and how to deal with 
that issue. This is likely due to the somewhat loose 
integration of computational modeling in their physics 
course. Presently, the length of time between exposures to 
VPython is too long, and students spend too much time 
relearning old programming skills. The course requires 
tighter integration of computational modeling into each 
assignment and modeling cycle. We have begun provid- 
ing scaffolded code and performing live coding exercises, 
both of which are best practices from computer science 
education. However, research in the field of computa- 
tional instruction is far from exhausted, and we have be- 
gun to develop our own studies of student thinking and 
practices. To a lesser extent, students have some trouble 
installing the software on their own computers; it requires 
several steps to install VPython and PhysUtil at present. 
We enlisted another team of Georgia Tech computer sci- 
ence majors to develop a web-based solution to streamline 
installation and updating the software (PhysKit) [15]. 

Resources for computational instruction are not 
widespread; most materials were developed by Georgia 
Tech's Physics Education Research group in conjunction 
with the classroom teacher. However, a virtual commu- 
nity of computational modeling instructors has begun 
building resources for math and science teachers inter- 
ested in introducing students to numerical computation. 
Many of these computational thinking resources are avail- 
able online. 4 Not all of these resources are tied to the 
Modeling Instruction curriculum nor are most resources 

4 http : //www. google . com/edu/ect/ 



physics related, but the support of such a community 
could produce additional high-quality resources. 

Students in these classes have had the support of a 
classroom teacher who works closely with University fac- 
ulty who are familiar with computational instruction in 
introductory physics. This support provides a strong re- 
source for the classroom teacher to collaborate on the de- 
velopment of materials, to identify and diagnose students' 
challenges, and to assess the implementation of the cur- 
riculum. Teachers who begin using these materials will 
not have the same close connection to the research group, 
but we will continue to provide support through sum- 
mer teacher development workshops on computational 
instruction at Georgia Tech. Moreover, a growing com- 
munity of users, developers, and researchers are start- 
ing to work with computational modeling in introductory 
physics courses. This community can provide additional 
support for early adopters, interested teachers, and, most 
importantly, our students. 
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